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Abstract 

We study the geometry of a codimension-one foliation with a time-dependent Riemannian 
metric. The work begins with formulae concerning deformations of geometric quantities as the 
Riemannian metric varies along the leaves of the foliation. Then the Extrinsic Geometric Flow 
depending on the second fundamental form of the foliation is introduced. Under suitable assump- 
tions, this evolution yields the second order parabolic PDEs, for which the existence/uniquenes 
and in some cases convergence of a solution are shown. Applications to the problem of prescrib- 
ing mean curvature function of a codimension-one foliation, and examples with harmonic and 
umbilical foliations (e.g., foliated surfaces) and with twisted product metrics are given. 
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1 Introduction 

A Geometric Flow (GF) is an evolution of a given geometric structure under a differential equation 
related to a functional on a manifold, usually associated with some curvature. The popular GFs in 
mathematics are the heat flow, the Ricci flow and the Mean Curvature flow. GFs play an essential 
role in many fields of mathematics and physics (see [1], [12], [H] for the heat flow and Ricci flow; 
and [5], P| - for foliated Riemannian manifolds). They all correspond to dynamical systems 
in the infinite-dimensional space of all appropriate geometric structures (e.g., metrics) on a given 
manifold. GF equations are quite difficult to solve in all generality, because of their nonlinearity. 
Although the short time existence of solutions is guaranteed by the parabolic or hyperbolic nature 
of the equations, their (long time) convergence to canonical geometric structures is analyzed under 
various conditions. A GF on a foliated manifold is called extrinsic, if the evolution depends on the 
second fundamental tensor of the foliation. 

One of the principal problems of extrinsic geometry of foliations reads as follows: Given a 
foliation T on a manifold M and an extrinsic geometric property (P), does there exist a Riemannian 
metric g on M such that T enjoys (P) w.r.t. g? 

Such problems (first posed by H. Gluck for geodesic foliations) were studied already in the 
1970's when D. Sullivan provided a topological condition (called topological tautness) for a foliation, 
equivalent to the existence of a Riemannian metric making all the leaves minimal. 

Let M""*"^ be a connected smooth manifold with a codimension-one foliation J^, and a vector 
field N transversal to J-. A Riemannian metric g on M is adapted if g{N, TT) = and g{N, N) = 1. 
Denote by M the space of adapted Riemannian metrics of finite volume for [M,T,N). 

The definition of the T-truncated {r,k)-tensor field S (where r = 0,1, and ^ is the TJ^- 
component) will be helpful, S{Xi, . . . ,Xk) = S{Xi, . . . ,Xk) (Xi £ TM). For example, g is the 
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J^-truncated metric tensor g on M, i.e., g{TJ^,TF) = g{TT,TT), and ^(A^, •) = 0. Another exam- 
ple is the second fundamental form b{X,Y) = g^V xY, ]^), where X,Y £ TT, and the Weingarten 
operator A{X) = -{VxNy of T w. r. t. g, extended on TM by b{N, •) = and A{N) = 0. Denote 
bm the symmetric (0, 2)-tensor dual to A"", hence, bm{X,Y) = g{A"'{X),Y). 

We call the Extrinsic Geometric Flow (EGF) on (M, J^, N) (see Section [3]) a family gt of adapted 
Riemannian metrics (i.e., gt{N,TT) = 0, gt{N,N) = 1) satisfying along 

dtgt = Sih). (1) 

Here S{b) is a symmetric J^-truncated (0, 2)-tensor expressed in terms of b. Recently, the author 
and P. Walczak studied EGFs on codimension-one foliations, see [ID] -[11] and [9l Parts 2 and 3], 

dtgt = Hbt), (2) 

where the symmetric J^-truncated (0, 2)-tensor h{b) is given by 

En— 1 
m=0 

and fm are either functions on M x M or symmetric functions of the principal curvatures of J-. 
The EGFs ([2]) reduce to the first order PDEs and yield quasi-linear hyperbolic PDEs. In the paper 
we introduce and study EGFs ([1]) of the form 

dtgt = V%h{bt). (3) 

The J^-truncated symmetric tensor h(b) can be expressed in terms of the first partial derivatives of 
the metric g. Therefore g i— V^hib) is a second-order differential operator. Hence, ([3]) reduce to 
the second order PDEs and yield quasi- linear parabolic PDEs (as will be shown in Sections [2HS]) • 
EGFs preserve the properties of codimension-one foliations: umbilical {A = Aid, see Proposi- 
tion [7]), totally geodesic (^ = 0), and Riemannian {V^X = 0). 

We show that EGFs serve as a tool for studying the following question on geometry of foliations: 
Under what conditions on {fm) o-nd [M^T^go) EGF metrics gt converge to one for which T 
enjoys a given extrinsic geometric property, e.g., is umbilical, totally geodesic, harmonic, etc? 

The structure of the paper is as follows. Section [2] collects main results of the paper: in 
Theorems [I]-[2] we study local existence and uniqueness of EGFs (introduced in Section [H); in 
Theorems [3l-|l] and Propositions [TH3] we apply the method to the problem of prescribing mean 
curvature of a foliation. Section [3 contains formulae for the deformation of extrinsic geometric 
quantities of a foliation as the Riemannian metric varies along J^. Section [8] contains examples for 
EGFs on foliated surfaces. Section [9] (Appendix) collects the necessary facts about parabolic PDEs 
and extend results of [9j about the generalized companion matrix. 

2 Main Results 

The power sums of the principal curvatures ki, . . . ,kn oi F (the eigenvalues of A) are given by 

r^ = k{ + ... + ki = tT{Ai) (j>0). 

The r's can be expressed using the elementary symmetric functions aj = Ylii< <i^«i " • • • " ki- 
where < j < n. Notice that ctq = 1 and tq = n. Set 7" = (ti, . . . , Tn) and = (ai, . . . , an). 

The choice of the right-hand side (in EGF equations) for h{b) seems to be natural; the powers 
bm are the only (0, 2)-tensors which can be obtained algebraically from the second fundamental 
form 6, while ri, . . . , r„ (or, equivalently, cJi, . . . , (t„) generate all the scalar invariants of extrinsic 
geometry. Powers bm with m > 1 in Q or ([5]) are meaningful; for example, EGFs 

dtgt = V'NTr{ht), (4) 

produced by the Newton transformation T^, see Example [H depend on all bm {fn < r). 
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Remark 1. Let Cn denote the Lie derivative w. r. t. N. Recall that C^f = N[f) for a function 
/ on M, and C^X = [N^ X] for a vector field. If a is any (0, j)-tensor, then 

Cn{Xi, ...,Xj) = N{a{Xi,. . . , Xj)) - ^ . a(Xi, . . . , X^.i, [N, X^], X^+i, . . . , X,) 

for any vector fields Xi, . . . ,Xj. The following EGF is similar to ([3]): 

dtgt = CNh{bt). (5) 
Since h{b) is a symmetric J^-truncated (0, 2)-tensor, £jv/i(&) preserves this property: 

CNh{b){N, X) = N{h{b){N, X)) - h{b){[N, N],X) - h{b){N, [N, Y]) = {X £ TM). 
Theorem [T] (and Theorem [5] in what follows) about EGF ([3]), are also valid for EGF ([S]). 
Define the quantities ^mij, depending on the principal curvatures {ki) of A, by 

= Y^lll k^kj-^-^ (1 < i < i < n, l<m<n). 
Remark that ^i-ij = 1 for all The ellipticity condition for EGFs takes the form 

En— 1 
, fm IJ-m;ij < 0, {I < i < j < n). (6) 

m=l 

In case of umbilical foliations (i.e., A = Aid), the condition ^ reads as ^ m fm\^~^ < 0. 

m=l 

Theorem [T] concerns the local existence and uniqueness of EGFs with fm of type (a) . 

Theorem 1. Let {M,go) be a closed Riemannian manifold with a codimension-one foliation J- , 
and let fm S C^iM x M). If (S^ is satisfied on M for t = 0, then there is a unique smooth solution 
9t of defined on a positive time interval [0,e). 

So, given fm G C^(M), if ([6]) holds for t = then there is a unique smooth local solution to 

En— I , 
^fm{t)VN bm- 

The generalized companion matrix Bn,i is defined in [9], see Section [9^21 Consider the matrices 

5n = 1 T /-(^) (^Mr-\ A={A,l A,, = ^Y.^~\^m-lfm,r,i^). (7) 

The next Theorem [2] about EGFs with fm of type (b) is based on Theorem [TJ 

Theorem 2 (Short time existence). Let {M^go) he a closed Riemannian manifold with a 
codimension-one foliation T and a unit normal N. Given fm £ C^(M"), suppose that at t = 
the condition ^ is satisfied and the matrix Bn + A, see is negative definite on M. Then there 
is a unique smooth solution gt of ^ defined on a positive t-interval [0,e). 

Example 1. (i) Let /(r^) > c> 0. Consider EGF 

dtgt = -f{'f^)V%bt. 

The corresponding system for r's (see Lemma[T]in Section [3] with S = f(T^)Vj\fb) is 

dtTi = \ fiT") N{N{ri)) +ai, i> 0, 
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where Oj = ^ N{f{'7^))N{Ti) — |/(t^) tr (VAr(^* ^)VAr(j4)). The above in-homogeneous parabohc 
system for r's along A^-curves (see Appendix), has a unique solution for < t < e. 

(ii) Consider EGF with fm of type (b), given along T by 

dtgt = V'^{f{l^)bt), 

where ^/(T^)id + A < {Aij = ^Tifj) for all t^((?) and a G M. The system for r's (see 
Proposition [9] in Section [9]) is parabolic quasi-linear 

dtn = - (^ /(T^) id + I(7^))iV(iV(r,)) + a^, i> 0, 

where a, = f /(7^) tr (V*v(^^-i)V*v(^)) - iV(/(7^)) iV(T,) - | r,E:,fe=i/,a6(^)iV(Ta)iV(Tb). The 
solution is unique and exists as long as ||7"|| is bounded. 

A foliation with vanishing ti is called harmonic; its leaves are minimal submanifolds. A foliation 
is taut if there is at least one metric on M for which is harmonic. The known proofs of the 
existence of such metrics use the Hahn-Banach Theorem and are not constructive. In particular, 
is taut if and only if there is a closed transversal curve that intersects each leaf see P], [7]. For 
example, a Reeb foliation on is not taut. 

Theorem [3] (based on Theorem [T]) helps us to find in some cases a one-parameter family of 
metrics converging to the metric for which is harmonic. 

Theorem 3. The following EGF, see on a closed Riemannian manifold {M^qq): 

dtgt = -2CNbt (8) 

admits a unique smooth global solution gt (t > 0). Moreover, if N -curves compose a fibration 
^ M ^ B, then the metrics gt approach to a smooth metric g as t ^ oo, and J- is g-harmonic. 

Remark 2. One may replace ([8]) in Theorem [3] by the PDE ()10p in what follows. By Example [7] 
(in Section [8|), the claim is false without condition 'W-curves compose a fibration". 

Let M = Ml X be the product of a Riemannian manifold {Mi,gi) and a circle {S^,g2)- For a 
positive function / G C^{M), the metric of a twisted product Mi x fS^ is given hy g = {f'^gi) © 92- 
The foliation F = Mi x {y} is umbilical with the mean curvature vector H = —V(log /)-*-. 

Proposition 1. Let Mi x jS^ be a twisted product of a closed Riemannian manifold {Mi,gi) and 
a circle {S^,g2). Then EGF (0) admits a unique smooth solution gt M for all i > 0, consisting 
of twisted product metrics on Mi x j^S^ . As t ^ oo, the metric gt converges to the metric g of the 
product (Mi,/i -^i) X (5\g2), where f{x) = J^i f{0,x,y)dyg. 

Example 2. Given f : T'^ ^ M+, let ^ f^^ be a twisted product of two circles {S^,gi). Then 
EGF ([8]) has a unique smooth solution gt £ M. for all t > 0, consisting of twisted product metrics 
on the torus T^. As t — )• cxd, the metric gt converges in to the flat metric g on T^. 

Proposition 2. Given f G C°°(R"), suppose that EGF on M, either (i) or (ii), 

(i) dtgt = -N{f{-^))gt, {ii) dtgt = -N{fi^))gt (9) 
satisfies the inequality for go (where functions Fk,'^^ o-fs defined in Proposition^: 

W El kFk-i{Ti)lk{T')>0, (ii) V" (n-A: + l)^fc_i(cTi)/,fc(^) >0. 
Then (0j admits a unique smooth solution gt defined on a positive t-interval [0,e). 
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Example 3. Notice that for / = § ri, EGF Q coincides with 



dtgt = --N{n)gt. (10) 
n 

We shah illustrate Proposition [2] for / = ^ (A; > 1) and f = ^ at {k > 1). Consider EGFs 

(i) dtg = --N{Tk)g, (ii) dtg = --N{ak)g. (11) 
n n 

If (i) Fk-i{Ti) > 0, (ii) ^'fc_i(cji) > 0, then the corresponding PDE is parabolic: 

(i) dm = N{Fk_,{n)N{n)), (ii) dtai = '1^^±1 N{^k-i{^i) N{ai)). 

Solving these systems and using relations Tfc = Ffc(ri) and ak = Fk{cri), see Proposition H] in 
Section [3l one may find a local solution of EGFs (jlip . 



Using Propositions [5] and [T2] (in Section I9.ip , one may show the following. 
Theorem 4. If for some real c > the following ellipticity condition holds: 

{i)c-'<y\ kFk^i{n)f,k{T') <c, {ii)c-^<y\ fn-/c + l)^fe_i(ai)/,fc(^) <c, 

then admits a unique smooth global solution gt (t > 0). Moreover, if N -curves compose a 
fibration ^ M ^ B, then gt approach to a smooth metric g as t ^ oo, and T is g-harmonic. 

A compact saturated domain D C M is called a {+)-fcd, if the transverse orientation of J- 
is outward everywhere on dD, and {—)-fcd if the transverse orientation of J- is inward on dD. 
P. Walczak and G. Oshikiri studied the question for codimension-one foliations (see [6], [15] etc): 

Which function F on M is admissible, i.e., can be the mean curvature of J- w. r. t. some 
Riemannian metric on M? It was proven that F is admissible if and only if 

F{x) > somewhere in any minimal [+)-fcd and F[y) < somewhere in any minimal {—)-fcd. 
In view of the integral formula by Reeb, see [9j, 

ridvol = 0, (12) 

M 

the condition J"^^ F dV = Q (for some oriented volume form dV on M) is necessary. 

In the next proposition we show how to produce in some cases a one-parameter family of metrics 
whose mean curvature (of J-) converges to a prescribed function. 

Proposition 3. Given F G C^(M) with the property Jj^j F dvolg = on a closed Riemannian 
manifold manifold {M,g), the PDE 

dtgt = --N{n-F)gt (13) 
n 

has a unique smooth solution gt {t > 0). Moreover, if each N -curve is dense in M, then lim r* = F. 



t—^ca 



3 J^-truncated variations of geometric quantities 

Let A^i C be the subspace of metrics of unit volume, and 

TT-.M^Mi, Trig) =9={ vol(M, g)-^/^ g) ® g' 
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the J^-conformal projection. Let gt € M (with < t < e) be a family of metrics with vol(M, go) = 1. 
Set St = dtgt- Observe that [X, Y]-^ does not depend on t. The volume form of gt evolves as, see 



^volt = ^(tr 5«)voU. (14) 

If St = St gt, where Sf : M — )■ R (i.e., gt are J^-conformal and gt is the D-truncated metric gt) then 
^ vol( = ^stvolt- Hence, the metrics gt = {4>t9t) ® with dilating factors cpt = vo\{M, gt)''^^"' , 
belong to A^i. 

Recall that the Levi-Civita connection V* of a Riemannian metric gt on M is given by 

2gt{V'xY,Z) = X{gt{Y,Z)) + Y{gt{X,Z))-Z{gt{X,Y)) 

+ gt{[X,Y],Z)-gt{[X,Z],Y)-gt{[Y,Z],X) (15) 

for all vector fields X, Y and Z on M. Since the difference of two connections is always a tensor, 
lit '■= 9tV* is a (l,2)-tensor field on (M,gt). Differentiation ()15p w.r.t. t yields the formula, 
see [m, 

gtiUtiX,Y),Z) = ^[iV'^St){Y,Z) + iV\.St){X,Z) - (V*^5t)(X, y)] (16) 

for aU X,Y,Z € r{TM). If the vector fields X = X{t), Y = Y{t) are t-dependent, then 

dtV'^Y = UtiX, Y) + V'AdtY) + V%^Y (17) 

Notice the symmetry Ilt{X, Y) = Ilt{Y, X) of the tensor 11^. 

Next, we find the variational formula for A, and apply it to symmetric functions tj and aj of A. 

Lemma 1 (see [9]). Let gt € Ai and S = dtgt ~ « J-' -truncated tensor. Then the Weingarten 
operator A of T (w.r.t. N) and the symmetric functions ti and Ui of A evolve by 

dtA = 1([A,5«]-V*v5«), (18) 

dtn = -^tr(^^-iV^5«), dtai = -UT{Ti_i{A)V%S^), i > 0. (19) 

For S = sg, where s S C^{M), we have 

dtA = -^N{s)id, dtTi = -^Ti^iN{s), dtai = -^{n-i + l)ai-iN{s), i > 0. 

Remark 3. Let (M, g = g (B g'^) he a Riemannian manifold with a codimension-one foliation J-" 
and a unit normal N. Given ip G C^(M), define a metric g = (e'^^g) © g'^. Then the second 
fundamental forms and Weingarten operators of w.r.t. g and g are related by, see [9], 

b = e'^'^{b-N{tf)g), A = A- N{^)(d. 

So, if N{ip) = (e.g., if ip is constant), then (i) b = e'^^ b; (ii) A = A; (iii) fj = tj. 

Proposition 4. Let gt £ Ai and dtgt = stgt- Then 

1 ,1. If] fk. 

rJ = F,(r,'), where F.(n) = rf + V 7I7X-7TT ' 3^ -f"'. (20) 



t , fx , , s (n — 1)! 1 I, v-^fe in — i)\ thj hi,-, 

.1 = where = ■ ^ E... (t - - .)! ' "f"' 

/or A; > 1. For small k, we have 

-F2(n) = - T"! + 4>2, F3{ti) = ^ + - (/)2 n + </)3, F4(ri) = ^ + (/)2 Ti + - (/"s n + (/>4, ete. 

^ , , n — 1 n -r / N (n— l)(n — 2) o n — 2 , 

2 n D n 

The t -independent functions 4>k = Ffc(O) and Vfc = ^fe(O) on M are determined att = by recursion. 
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Proof, (i) We will prove ()20p by induction. Notice that the last term (with i = k) m ()20p is (pk- 
By Lemma [T| the system for r' reads as 

dtTi = -'^N{s), dtT2 = -nN{s), dtTk = -^Tk-iN{s), ... (22) 

By (l22])i we have 9t(rf) = -nriiV(s). Substituting N{s) into (l22])2, we find 5t(r2 - ^rf) = 0. 
Hence, T2 — -^rf = 4>2 is t- independent, and F2(ri) = ;^ + 02, see (PU|) for k = 2. 
Assume that (|20p holds for some k > 2. Using ([22]) 3, we find 



k + 1 ^ , , , k+1 



1 . 1 k\ 



= -— F.(n) iV(») = -— ri + . A rf-'J iV(<,). (23) 

By (f22]) i we have 5t(r^'''"'^) = _ "(■^^-'-^ j-^ ^{s). Substituting this into ([23]) . we obtain dtTk+i = 
and hence (pO]) for /c + 1. 
(ii) We will prove ([2T]l by induction. Notice that the last term (with i = A;) in ([2T|l is V'fc- 
By Lemma [U the system for cr' reads as 

atCTi = -^iV(s), 9^(72 = -^^f7iiV(s), dtak = -'^^^^^Ok-iN{s), ... (24) 

By we have 5t((Tf) = -nai N{s). Substituting N{s) into (^2, we find 5t(c72 - ^ (^1) = 0- 

Hence, ct2 — af = tjj2 is t-independent, and ^'2(0'!) = ^1 ~^ ^2, see (^T]) for k = 2. 

I 



Assume that (|2ip holds for some k>2. Using ()24p ci. we find 



n — k 



(n-1)! 1 ^ 1 (n-i) 



.fc!(n-fc)! n^-T ^1 + Z^.=i " (A; - i)!(n - fe)! ""i"'] ^^'^^ ^^^^ 

By (p^ i we have dt{cri'^^) = _ "(■^+^) q-j N{s). Substituting this into ([25]). we obtain dtat+i = 
9t*fc+i, and hence ^ for A; + 1. □ 

Remark 4. We may set Fq = n, Fi{ti) = ti and (pi = 0. Differentiating the relation = -Ffc(ri) 
with A; > 1, see Proposition HI and replacing dfTi = — ^ ^(s) by Lemmadl we find 

dtn = F^{n) dtn = FUri)Nis). (26) 

Comparing ([26|1 with dfTk = — | t^^i N{s) = — | Ffc_i(Ti) N{s), see LemmalU we conclude that 

F^(ri) = -Ffc_i(Ti) ^ Fk{Ti) = -['Fk-iix)dx + cj)k. 
n n Jq 

Similarly, for functions ^'^(cri) in ([2T]) we set ^'o = 1, ^i(cri) = o"! and V'l = 0, and obtain 

^fc(f^l) = ^fc-l(fTl) ^ ^fc(fTl) = / ^'fe_i(x)dx + Vfc- 

71 n Jo 

Lemma 2 (see [9]). The vector field Z = V^j^N is evolved by gt € Ai with S = dtgt as 

(?) dtZ = -S^{Z), (ii) dtZ = -sZ for S = sg. (27) 

In particular, all variations gt G Ai preserve Riemannian foliations. 

We will use the following condition for convergence of evolving conformal metrics (for a general 
formulation see [H Appendix A] ) . 

Proposition 5. Let dtgt = st gt (t > 0) be a one-parameter family of J- -conformal Riemannian 
metrics on a closed manifold M with a codimension-one foliation T and a unit normal N. Define 
a function v{t) = supj^,^ I'^tlgC*) '^'^'^ assume that v{t)dt < 00. Then, as t ^ 00, the metrics gt 
converge in C° to a smooth limit Riemannian metric goo ■ 

Proof. Our assumptions ensure that the metrics gt converge in to a symmetric (0, 2)-tensor 
goQ. Notice that the metrics are uniformly equivalent: c^^go < gt < ego for some c > and all 
t > 0. Hence, goo is positive definite. □ 



7 



4 The Extrinsic Geometric Flow 

We study two types of evolution of Riemannian metrics on (M, T), depending on functions fm {m < 
n), at least one of them is not identically zero: (a) fm G C^(M x M), (b) fm = /m(T^) G C^(M"). 

Definition 1. Given functions fm of type either (a) or (b), a family gt, t £ [0,e), of adapted 
Riemannian metrics on {M,T,N) satisfying along T the PDE see also ([5]), will be called an 
Extrinsic Geometric Flow (EGF). The corresponding normalized EGF is defined along T by 

dtgt = VNKh)--r{t)gt, where r{t) = [ N {tr h{b)) d volt / yoI{M, gt). (28) 

Lemma 3. Let gt be the solution to EGF Then the Weingarten operator A of T w. r. t. gt 
satisfies 

= Er=o il ^^(^")] - N{N{fm))A^ - N{fm)V%iA"') - fmV%V%{An) . (29) 

Proof. Substituting = V^/i(^) into ((IHD, see ([3]), we obtain i^. □ 
Let gt be a family of Riemannian metrics of finite volume on {M,F). Metrics gt = {(ktgt) © gt 
with (f)t = vol{M, gt)~'^^"' have unit volume: Jj^jdvolt = 1. The next proposition shows that 
unnormalized and normalized EGFs differ only by rescaling along T . 

Proposition 6. Let {M,T) be a foliation, and gt a solution (of finite volume) to EGF with 
^(^) — Y^'fn=o fmCT^) bm- Then the metrics 

gt = {4>tgt)®g^, where (t)t = vol(M,54)-2/", 

evolve according to the normalized EGF 

dtgt = y%h{bt)-{pt/n)lt, where pt= [ N {tr h{b)) d volt / vol{M,gt). (30) 

Proof. By Remark El fj = tj and h(A) = h{A) (for metrics gt and gt, respectively). Hence, 

tr h{A) = tr h{A) = V""^ fm{^) Tm- 

^ — ^m=0 

From (HID with S = V%h{b) we get the derivative of the volume function 

d d C 1 /* 

— vol{M,gt) = — dvolt = - N{tY h{A))dvolt. 
at at Jm i 

Thus (t>t = vol(M,5t)-2/" is a smooth function. By Remark [3l we have h{b) = (pt • h{b). By the 
above and gt = 4>t9t, we have 

dtgt = (t^tdtgt + 4>'t9t = h(b) + (l}'t/(l}tgt- 

Using (fH|) and dvolt = (j)^^"^ dvolt, we obtain 

dtvolt = dtict^l volt) = - (j)^ <P't volt +- N{tv h{A)) volt = - g + iV(tr h{A))j vol*. 

Let pt be the average of A^(tr h{A)), see (pOj) . From the above we get 

= 2^ / d^lt = [ (n^ + N{tr h{A)))d^lt = n^ + pt. 
dt Jm Jm ^ n ) <Pt 

This shows that ptjn = —(^[/(pf Hence, g evolves according to (PU)) . □ 
Obviously, EGFs preserve Riemannian foliations. Totally geodesic foliations are the fixed points 
of ([3]) . We shall show that EGFs preserve the umbilicity of J^. Define the function of one variable by 

^ — ^m=0 
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Proposition 7. Given go on (M,J^), let (?t (0 < t < e) be a unique smooth solution to the EGF 
with fm of type (b). If J- is umbilical for go and > then T is umbilical for any gt. 

Proof. Since is go-umbilical, we have = '^o id for some function Aq : M — t- R. Hence, 
/i(Ao) = — V'(Ao) id. Assume that J- is umbilical for all gt, then 

^t = Atid, h{At) = -ip{Xt)(d, 

where Aj : M — t- R are smooth functions. In this case, by Lemma [1] with s = — '0(Aj), we have 

dtXt = \ N{N{i;iXt))) = \ iV(Vy(At)iV(Ai)). (31) 

Since ^^'{X) > 0, the PDE (f3T]) is parabolic and admits a unique local solution, see Appendix. Now, 
([3]) reads as 

dtgt = -N{i;{Xt))lt, (32) 

from which we find gt = exp(— J^^^ Nlip^Xs)) ds)go along J^. The Weingarten operator of certain 
metrics gt is conformal: At = A^ id. Thus, gt is the solution of ([3]). Our assumptions ensure that a 
solution is unique, hence gt = gt for alH. □ 

5 Searching for power sums 

Although EGF ([3]) consists of (second order) non-linear PDEs, the corresponding power sums 
Tj (i > 0) satisfy an infinite quasilinear system. 

Proposition 8. Power sums {rjjjgN of EGF (0) with fm of type (a) satisfy the infinite system 



dtn = E " , N{N{n^^^^)) + a„ (33) 

2 ^-^m=i 1 + m — 1 



n-l 

where ai = -^ [N{N{fm))Ti+m-i- fmtr iV%iA''-'^)V%{A''"^))]. The n -truncated system ^ is 

m=l 

dt^ = - E""' ^ fm iBn,ir-' N{Ni^)) + a, (34) 

where a = a(T^, A^(t^)) is a smooth vector-function. 

Proof. Since we substitute 5 = V^j^h{b) into ([121) i ^i^d obtain 

dtTi = -'-tr (A*-iv5vV^( V""' /^^"^)), i > 0. (35) 

2 ^ — ^m=0 

The desired system ()33p follows from the above and the identity 

'^"""^ tr {A'''V%A'^) = tr {V%A^+"'-') = iV(r,+„-i). (36) 
m 

By Proposition [13] (see Appendix), the system ([33]) is equivalent to ([M]) . □ 
Proposition 9. Power sums {TjjjgN of EGF with fm of type (b) satisfy the infinite system 

i r ■sr-^n-l mfm iKT/T^T/ w 

i + m 

where N{fm) = J2s fm,sN{Ts) and 

. n—l „ n 
«^ = -'l^y,\--^N{fm)N{Ti+m-l)-fmtT (V5v(^*-')V,v(^™)) +r.+^-l V fm,abN{ra)N{n) 



= -T^^ ^ iV(iV(r.+^-i))+E""n^^+>n-iE" 1 (37) 

2 L ^ — ^m=l I 4- m — L ^ — ^m=0 ^ — ^s=l J 



m=l a, 6=1 
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The n-truncated system ( fgTj j is 

dtl^ = -{Bn + A)N{N{l^)) + a, (38) 
where a = aij"^ , N {1^ )) is a smooth vector-function, the matrices and A are defined in ([^. 

Proof. It is similar to proof of Proposition [H Substituting S = V]\fh{b) into (fT9]) i . see ([3|), we 
obtain ()35p . The desired system (|37p follows from the above and ()36p . By Proposition 1141 the 
system (j37|) is equivalent to ([38]) . □ 

Problem. Express Oj (i > 1) in Propositions [8] and [9] through the scalar invariants of A and V^A. 

Example 4. Among EGFs (l3|) with fm of type (b), the Newton transformation flow, see (jH), 
corresponds to /m = (— l)'"crr_m- Newton transformations Tr{A) are defined as 

rr(^) = y"'^ (-l)V^_iA^ = arid-ar_i^ + ... + (-l)M'^ (0 < r < n), (39) 
see ^J. Since tr Tr{A) = (n — r) ar, the corresponding to dH) normalized flow is 

dtgt = V*;v rr(6i) - - r(t) 5i , where r{t) = {n - r) [ (a,) d voU / vol(M, gt). (40) 
For EGF denoting Tj = Tj(A), by Lemma [H we have 

dta, = -^[N{triT,^,.V%Tr))-tr{iV%T,^i){V%Tr))] = ^ tr ((V^ T,_i)(V^ T,)) 

N{ar^j tr (Ti^i • V*;v (^')) + A^(^Tr-j) tr (r,_i • ^^■)) . 
One may express the above (by induction) through it's and their A^-derivative. 



(-1) 



6 Local existence and uniqueness of metrics 

We shall prove our results concerning the existence/uniqueness of EGFs ^ and their corollaries. 
For brevity we shall omit the index t for time-dependent tensors A, a, bj and functions Ti,ai. 

Generally, EGFs with fm of type (a) are solvable under additional conditions (e.g., for A and 
the auxiliary functions fm)- The following theorem concerns EGFs with fm of type (a) and is 
essential in the proof of Theorem [2] about EGFs with fm of type (b). 

The coordinate system described in the following lemma, see |2i Section 5.1], is here called a 
biregular foliated chart. 

Lemma 4. Let M be a differentiable manifold with a codimension- 1 foliation T and a vector field 
N transversal to T. Then for any q ^ M there exists a coordinate system (xq, xi, . . . x^) on a 
neighborhood Uq <Z M (centered at q) such that the leaves on Uq are given by {xq = c} (hence the 
coordinate vector fields di = (9^., i > 1, are tangent to leaves), and N is directed along do = 
(one may assume N = do at q). 

If (M, J^, g) is a foliated Riemannian manifold and A^ is the unit normal to J-, then in biregular 
foliated coordinates {xo,xi, . . . , a;„), one may present the metric g in the form 

g = goo dxl + ^ . gij dxidxj , (41) 

where one may assume goo = along the axis {xo,0, ... ,0). Let g^^ be the entries of the matrix 
inverse to (gij) and gij^k the derivative of g^j in the direction of dk- We have bij = g{VQ^dj, N). 
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Lemma 5 (see [9]). For a metric (^Tjj in biregular foliated coordinates (of a codimension-one 
foliation T) on {M,g), one has N = Oq / y/goo (the unit normal) and 

r|o = (1/2) J]^fc,o/^ no = -(l/2) j;^5oo,s/\ ro. = -5,,-o/(2<7oo), 

bij = r^ -v/goo = — dijfl I \/goo (the second fundamental form), 

A{ = -r^n/Vgoo= o .j— V] gis.og''-^ (the Weingarten operator), 

{bm)ij = A^... Al2~^A^"^gjs, (the m-th "power" of hj). 



52 

-I \i 



Sir. 



n = (75-7==) V, , , , 5ri.2,o---5n_is„o/""' •••5 
In particular, n = Yl 9rs,o 

^ r, s 

Proof of Theorem [TJ Given g G M, by Lemma [H there exist biregular foliated coordinates 
{xo,xi, . . . Xn) on Ug C M (with center at q) and the metric has the form (jH]). Then the unit 

— 1/2 

normal to T on Uq is N = g^^ do. Denote Fij := h{b)(t, x, gab,'ipab)ij and V'ab = d^gab- The EGF 
system ([3|) with of type (a) has the following form along A^-curve 7 : x 1— )• 7(x) (on Ug): 

dt9ij = 9oo^'^ {doFij - LQjFfcj - TQjFik) 

= g-^'^ {doFij - i ^^^^{i^isg"'Fkj + ^isg'^Fik)) {i,j = l,...,n), (42) 

where ^oo = 1 along the axis 70 = (a^o, 0, . . . , 0). Certainly, EGF system ^ has a shorter form 

— 1/2 

dtgij = 9oo doFij{t,x,gab,Tpab)- (43) 
For example, if = for m > 2, then Fij = foQij + fibij {i,j > 0), and ()32]l reads as: 

dt9ij = -^9oo^^fi dio 9ij + aij{t, x, gab,'^Pab), 

which is parabolic when /i < 0, see also ([6|). Using Theorem A completes the proof in this case. 
Now let /m 7^ for some m > 2 (e.g., general /m). Computing the derivative 

9o^ij = ^009ab + i'ab + 

^a<6 dlpab ^a<b dgab Oxq 

we rewrite (j42p as the quasi-linear system 



k9ij = 9oo yZ <^h WT^ ^lo9ab + Cij{t, X, gab,ipab)- (44) 



In general, the following square matrix of order ^ n{n + 1) is not symmetric: 

{Qp. . -| 
jr^\ {i<j, a<b). 

We claim that ()44p is strong parabolic. If we change the local coordinate system on M, the 
components Fij (i < j) and V'afe (o < &) at g will be transformed by the same tensor low. Hence, 
d^F is a (1, l)-tensor on the vector bundle of symmetric (0, 2)-tensors on TT. So, d^F{q) can be 
seen as the linear endomorphism of the space of symmetric (0, 2)-tensors on TgT. The order of 
indices in matrix representation of d^F is [1, 1], [1, 2], ... , [1, n], [2, 2], [2, 3], . . . , [2, n], . . . , [n, n]. 
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The strong parabolicity is a pointwise property, so can be considered at any point q € M in a 
special biregular foliated chart around q, for example, such that gij = 6ij and Aij = ki5ij at q for 
t = 0. (Recall that (ki) are the principal curvatures of J-). In this chart, our calculations show 
that the matrix d^F is diagonal, so has real eigenvalues at q, and 



dFi 



dip, 



ab 



E 



m>l 



/m(g,0) 



By ([6]), the system ()44p is strong parabolic: the constant in the condition ()55p . see Appendix, is 

c = infinf|> fm{q,^) ^J■m■,^j\ > ^■ 

q i,j ^ ^ — ^m=l 



By Theorem A (in Appendix), given q £ M there exists a unique solution to (|44p which is defined 
in Uq along the A^-curve through q for some time interval [0, e^) and satisfies the initial conditions 
gij{0,x) = go,ij- By Theorem A, the value Eg continuously depends on g G M. The claim follows 
from the above and compactness of M. □ 

Example 5 (see |9]). For h{b) = 62, the matrix d^F[q) in an orthonormal frame at any point is 



'2)ij 



dip, 



ab 





2Vl2 


2tPl3 








- 


■012 




l/'23 


■012 


013 









V'll+'033 





V'12 


V'13 





2Vl2 





2 V'22 


2t/'23 








i>13 


■012 


1p23 


■022+ -033 


V'23 








2 013 





2 023 


2^33 - 



with the order of indices [1, 1], [1, 2], [1, 3], [2, 2], [2, 3], [3, 3]. At g for t = (i.e., ipab = 0, a^b and 



Ipa 



ka) it is diagonal with the elements //2;afe = ka + kb- For instance, let Aq = [i, 1, 1, ^, 1, h] be 



the diagonal matrix of order six. Then the matrix Ai = AQd^{b2) is symmetric, and the system 
for h{b) = 62 is "symmetrizable" : A^dtgij = Aid^^gij + {low-order terms}. 

Remark that for EGFs with general fm of type (b), the solution procedure requires additional 
assumptions (e.g., for the auxiliary functions fm and the matrix Bn)- 

Proof of Theorem [2J Let Aq and 7" be the values of extended Weingarten operator and power 
sums of the principal curvatures ki of (the leaves of) F determined on (M, F) by a given metric 50- 

a) Uniqueness. Let g^'^ and gf''^ be two solutions to ([3]) with the same initial metric g^. Functions 
T^*'^ and 7"*'^, corresponding to g^^'^ and g't'\ satisfy (j37p and have the same initial value r^''. 
By Theorem A (see Appendix), a solution of the parabolic system ()37p is unique, and we have 
T^*'^ = T^*'2 = T^* on some positive time interval [0,ei). Hence, both g^^^ and gf'^ satisfy EGF 
equation ([3]) with of type (a) with known coefficients fj{p,t) := /j(7"*(p)). By Theorem [U 
gl^^ = g^^ on some positive time interval [0, £2)- 

b) Existence. By Proposition [U] and Theorem A, (j37p admits a unique solution 7"* on a positive 
time interval [0, ei). By Theorem [H Q with fm of type (a) with known functions fj{-,t) := 
/j(T^*, t) has a unique solution g^ {gQ = go) for < t < e*. The Weingarten operator A^ {Aq = Aq) 
of {M,F,g^) satisfies (f29|) . hence the power sums of its eigenvalues, 7^*'* (1^^'* = 7^''), satisfy 
(|33p with fm of type (a) with the same coefficient functions fj. By Proposition [8] and Theorem A, 
the solution of the problem is unique, hence r^* = r^*'*, i.e., r^* are power sums of eigenvalues of 
A'l . Finally, g^ is a solution to ([3]) such that 7" * are power sums of the principal curvatures of F 
in this metric. □ 
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7 Proof of Theorem [3] and Propositions [TH3] 

Proof of Theorem [3l In biregular foliated coordinates, ([8]) has the form of parabohc PDEs 




see also ()33]l . One may assume goo = 1 along the central A^-curve 70 (s), hence dtQij = 9oo 5ij 
along 70- This linear system has a unique global solution gij{t) for all f > 0, which converges to 
g = (jjij) as t — )• 00 (along any A^-curve), see Appendix. Hence there is a global solution gt {t > 0) 
on M. Since A^-curves compose a S'^-fibration, from dt{detgij) = ^{tv gS) det gij for S = —2£]\fb, 
see (|14p . we conclude that g is positive definite and the limit g = lim gt is a smooth metric on M. 

The corresponding system for r's, see ([3^ with i?^ = ^ id, is diagonal parabolic, 

dtn = N{N{Ti))-itv {Vn{A'~')Vn{A)). 

In particular, we have dtTi = N(N{ti)); from this easily follows that N(ti) satisfies the heat 
equation dt{N{Ti)) = N{N{N{ti))). Hence, the limit N{fi) = const along A'"-curves. Since fi is 
bounded on a compact manifold M, we obtain N{fi) = 0. By the integral formula Jj^i^iTi) — 
ff] dvol = 0, see [9], we conclude that ri = on M. □ 

Proof of Proposition [TJ Define (/) : M — )• M by the equality f = . We have ri = —N{cl)). 
The evolution dtg{t) = {s{t) gi{t)) © preserves twisted product structure. Denoting gt = gi{t), 
we obtain g{t) = (/(t)^5i) © 92, where 52 does not depend on t. Assuming f{t) = e^^^\ we find 
= e2<^W5o- In this case, dtg{t) = 2e2<^W9t(/>(i)5o- Hence, dt(l){t) = i s{t). For s = -^N{ti), EGF 
equation ([8]) reads as 

= {l/n)N{N {(/))). (45) 

Replacing t — )• t/n, one may reduce ([1^ to the heat equation. Hence, see Proposition!^ there is a 
unique solution cf){t,x,y) for all t > 0, satisfying lim (f){t,x,y) = ct){x). Certainly, / = e"^. □ 

Proof of Proposition [2l By Remark [H we have FI{ti) = -F]._i{ti). By Proposition H] with 
s = -NifC^)), we find 

By this and Lemma [U we conclude that ti satisfies the quasi- linear PDE 

dtn = N{a{T')N{n)), where a = ^ ^^^^ A; /,fc(7^)Ffc_i(ri). (46) 

This PDE is parabolic when a > (see Section [9. ip that proves (i). The proof of (ii) is similar. □ 

Proof of Proposition [3j By Lemma [1] with S = s g and s = —N{ti — F), we conclude that the 
functions ti — F and N{ti — F) satisfy the heat equation (see Appendix) along A^-curves, 

n 2 

dt{n-F) = dtTi = --N{—N{n-F)) = N{N{n-F)), 

2 n 

dt{N{n - F)) = dt{N{n)) = N{dtn) = N{N{N{n - F))). 

There is a unique solution, r{, for all t > 0; moreover, r{ — >• fi uniformly as t — !• 00, with the 
property N{fi — F) is constant along A^-curves. Since M is compact, we find that N{fi — F) = 0, 
hence fi — F is constant along A^-curves. With known r^, (jlSp also has a unique global solution 
gt = go exp(— /q N{tI — F) ds) for all t > 0. Let us apply the normalized EGF 

dt~gt = -[-N{n - F) + -r{t)]lt, r{t) = -2 [ N{n - F) d^U / vo\{M , ~gt) ■ 
n n Jm 
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By Remark [31 the mean curvature is the same for both EGFs: r* = f^. Let us differentiate the 
functional Jpit) = /^^ Fdvolt, using with st = -| N{ti - F) - f r{t): 

Mt) = - [ N{n -F)d^lt-l r{t)Yo\{M,~gt) = 0. 
Jm ^ 

By conditions, Ji?(0) = 0. Hence Jpii) = 0, i.e., F has zero average on (M, ^t). Since A^-curves 
are dense in M, we have f\ — F = C for some C G M. Thus, r * — F = C + o(l/t). Integrating over 
(M, gi^), we obtain 

[C + o(l/t)]vol(M,gt) = / {T\-F)d^\t= [ rld^h. 

Jm jm 

Since Jj^^ t\ d volt = (the ^j-average of r* is zero), see (fT2]) . and vol(M, ^rt) = 1, we have C = 0. □ 

Proof of Theorem [4l By Proposition [21 ri satisfies parabolic PDE (jl6l) . By Proposition [121 there 
is a global solution t\ [t > 0) and we have r* = Fj(T*) for i > 1. The metrics gt {t > 0) are given 
by the formula gt = goe~ f^o fi"^ "^'^'^ . Moreover, there is a limit (as t — )■ oo) which is constant 
along A^-curves. By integral formula J^[N{ti) — (ri)^]dvol = 0, see [9], applying for gt, we have 
rf° = on M. The rate of converging — )• -Fj(O) (i > 0) is exponential, see Proposition \T1[ 
By Proposition [5l see also the proof of Theorem [3l there is a limit metric ^oo = lim ft ■ The case 
(ii) is similar. □ 

8 EGF on foliated surfaces 

In this section, (M^, (70) is a surface equipped with a transversally orientable foliation by curves. 
Given ip G C"^ (M) , let a 1-parameter family of Riemannian metrics ft (0 < t < e) satisfies the PDE 

dtgt = -N{tP{Xt))gt (47) 

see ([lOD and ([32]). Applying 1^ with 5 = -N{ip{Xt)) gt, we have ^ (rfvolt) = A^(?/'(Aj)) dvol*. 
Hence, the volume of a closed surface {M'^,gt) satisfies the equation 

d If n f 

- vol{M,gt) = -- Nit^{Xt))dvo\t = -- A* V(At) dvoU . (48) 



Proposition 10. IfO< Co(A) < '(/''(A) < C'i(A) < 00, i/ien on a surface {M'^,F) has a unique 
global smooth solution gt (t >0). 

Proof. The (^j-geodesic curvature At of F (w. r. t. A^ - the gt-unit normal to F) satisfies the PDE 

dtX = ^Nii;'{X)NiX)), (49) 

see (|3ip . By Proposition [12] in Section [9. H ()49p (with Aq determined by go) admits a unique smooth 
solution Xt for all t > 0. In this case, the EGF ()47p has a unique smooth solution gt on a surface 
for all t > 0, moreover, gt = go exp(— N{ip{Xt)) dt). □ 

Example 6. Let ip = 2X, i.e., dtg = —N{2Xt) gt, and conditions of Proposition [10] are satisfied. 
Then ()49p is the heat equation dtXt = N(N(Xt)) along A^-curves, whose solution is known. The 
EGF solution is gt = 5oexp ( — A^(Ag)d^) {t > 0). For closed manifold M, by ([ISp . we have 
^vol{M,gt) = /o A| < 0, hence the function vol{M,gt) is decreasing. 

Proposition 11. The Gaussian curvature of the EGF ^7\ ) on {M'^,F) is given by formula 

Kt = div (exp (^* Ar(V'(A5))d^)V^Af) + N{Xt) - Xl (50) 
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Proof. Notice that tj = nX'. By Lemma [2] with S = —N{ip(\t))gt (and n = 1), we have 
dtiVj^N) = N{^p{Xt))V'j^N. Integrating the above ODE yields 

V%N = ei,p(^I^N{^P{X^))dc)v%N. (51) 

In our case here, K = Ric(A^, A^), therefore (|5U|) is a consequence of (|51|) and the identity for 
codimension-one foliations (see |9l Proposition 1.4] with r = 0), 

div(V7viV + TiN) = Ric(iV, N) + T2- {nf. 

Indeed, one should apply the identities divA^ = — ri and div(TiA^) = N(ti) — (ti)^. □ 

Example 7. Let a function / G C^(— 1, 1) has vertical asymptotes x = ±1 and exactly one strong 
minimum at x = 0. For example, / = jQ{e^ /(i-^ ) _ i^. Consider a Reeb foliation J- in the strip 
n = [—1, 1] X M (equipped with the flat metric go) whose leaves are 

L± = {x = ±1}, Ls{x) = {{x, f{x) + s), |x| < 1}, where s E M. 

Identifying the points (— l,y) ~ (1)1/)) we obtain a foliated cylindrical surface x M. Next, 
identifying the points (l,y) — (l,y + n) for n G Z, we obtain a foliated flat torus 11' = 5^ x S^. 

Denote a{x) the angle between the leaves Lg and the x-axis at the intersection points. That is, 
/ and a are related by 

/'(x) = tana(x) ^ cos a = [1 + (/')^]~^/^ sina = /' [1 + (/')^]~^/^. 
The geodesic curvature of Ls is 

Ao(x) = /"(x)[l + (/'(x))2]-3/2 = a'{x) ■ I cosa(x)|, |x| < 1. 

The tangent and unit normal to J- vector fields (on the strip) are X = [cosq(x), sinQ(x)], N = 
[— sina(x), cosa(x)]. The unit normal N at the origin is directed along y-axis. We use X and N 
to represent the standard frame ei = cosa(x)X — sina(x)A^, 62 = sina(x)X + cosa(x)A^ in 11. 

Taking a{x) = ^x, we obtain the Reeb foliation with Ao(x) = ^ cos(^ x) > 0. 

Let the Riemannian metrics gt for t G [0,e) satisfy ()47p on 11'. Hence gt = ^oe~^*, where 

Ut{x)= f N{i^{X^{x)))d^ = -sma{x) f d.,iP{\^{x)) d^. 
Jo Jo 

So, gt{X,X) = e~^\ gt{X,N) = 0, and gt{N,N) = 1. The scalar products gii{t) = gt{ei,ei), 
512 (i) = 54(61,62) and 522 (*) = 54(62)62) are 

511 (0 = sin^ Q + cos^ ae~^\ guit) = sin a cos a(e~^* — 1), 522 (*) = cos^ a + sin^ ae~^K 

In particular, \gt\ = gii{t)g22{t) — 5i2(i)^ = 6~^*. The A^-curves 7(5) = (x(s),y(s)) satisfy the 
system of ODEs ^x = — sina(x), -^y = cosa(x). From the first of ODEs we deduce the implicit 
formula s = — J^"^^^ sin'a(g) local diffeomorphisms ips {s > 0) of the interval (—1, 1). (Indeed, 
moving along A^-curve by the length s from a point (x, y) to a point (x', y'), we find that x' = (^^(x)). 

Let tp{X) = 2A, see Example [H By Proposition [TOl the EGF dtgt = —2N{Xt)gt admits on 11' a 
unique global solution gt (t > 0), and (by Example [6]) the function vol (11', (7^) is decreasing. Hence 
5t = 5o6~^% where Ut = — sina(x) dxX^{x) d^, and A evolves along A^-curves with arc-length 
parameter s G M by the heat equation dtX = dg^X. In terms of x, we have dtX = sma{xyd'^xX. 
This heat conduction PDE in a material with variable "thermal conductivity" sina(x)'^ on a circle 
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is parabolic everywhere except one (singular) point a; = 0. Using Example [8] (i), we write solution 
along A^-curves as 

Xt{Ms)) = I^HMO)G{t,s,Od^, where G(t, s, = e"(^-«)'/(^*). 

Since Ao((/>o(-s)) quickly as |s| — )• oo, we have uniform convergence Xt{x) — )■ as i — )• oo for any 
X £ (—1, !)• Along the leaves L±oo the metric does not change, hence Ai(ibl) = for all t >0. 

It is known that a curve on a surface with A = is a geodesic. Prom the above we conclude 
that gt do not converge as t — )• oo to a Riemannian metric on a torus H', see Remark [2] (otherwise 
T becomes a geodesic foliation, that contradicts to the fact that Reeb foliation is not geodesible). 
Since gt do not depend on the y-coordinate, the Gaussian curvature is given by (see ^) 



\9t\ ""^y/lgtl 

Substituting 522 (*) = cos^ a + sin^ a ■ e~^* into the above formula yields 



2e-^'Kt = -(2cos(2Q)(a)^ + sin(2Q)a")(e"^'-l) + sin2ae-^'5^^f7f 

+ ^ sin(2a) a'e-^'d^Ut - ^ sin(2a) a'd^Ut - ^ sin^ a e'^'id^Utf. (52) 

Because a(0) = and C/t(0) = 0, by ([52]) one has i^t(O) = 0. For a = ^x and small x we have 

e-^*irt = ^7r3l4(0)x + o(x), 

o 

where 14(a;) = /q 92;A^(x) d^. Thus Kt{x) (for t > 0) changes its sign when we cross the line x = 0. 

9 Appendix: PDEs 

In Section 19.11 follow [12] we discuss parabolic PDEs and prove Proposition [12] about the quasi- 
linear heat equation. In Section |9.2] using the generalized companion matrix, defined in [9j, we 
prove Proposition [TJ] about infinite system of parabolic PDEs. 

9.1 Parabolic PDEs in one space variable 

Now, we recall some facts about parabolic systems of quasi-linear PDEs. 

Let A = {aij{t,x,u)) be an n x n-matrix, a = {ai{t,x,u,dxu)) - an n-vector, and t,x EM. 
Consider the quasilinear system of PDEs, n equations in n unknown functions u = (ui, . . . , Un) 

dtu = yl(t, X, u) d^xU + a{t, x, n, dxu). (53) 

PDEs ([53p are homogeneous if a = 0. When A depends on t and x only, the system is semilinear. 
When A and a are functions of t and x, the system is linear. The initial value problem for (|53p 
with given smooth data A, a and uq, 

n(0,x) = no(x), (54) 



consists in finding smooth function u{t, x) satisfying ([53p - ([54 

The parabolicity condition for ([53]) says that there is c = const > such that 



{Av,v) >c{v,v), VuGM". (55) 

Theorem A [12, Ch. 15, Proposition 8.3]. Suppose that A of class C°° satisfies the parabolicity 
condition Ii55\) . Then the system i53 \) - lji54\ ) with uq € H^{M) (the Sobolev space) for some s > 1, 
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admits a unique solution u € C([0, T), H^{]R.)) n C°°((0,T] x M), which persists as long as \\u\\c^ 
is bounded, given r > 0. 

One may apply this and the maximum principle to obtain the global existence and uniqueness 
result for a scalar parabolic equation (i.e., n = 1 and An = k{u) > 0) on S^. 

The following proposition is standard, for the convenience of the reader, we give its proof. 

Proposition 12. Let a function k{u) G C°°(M) (the thermal dijjusivity) satisfies 

< ci < k{u) < C2 < OO 
for some real c\ < C2- Then the quasi-linear heat equation on a unit circle 

dtu = d^{k{u) d^u), n(0, x) = uq{x) e C^{S^) (56) 
admits a unique solution u G C°°([0, oo) x S^). Moreover, there exists lim u[t,x) = Uoo £ and 

t~>-oo 

for some real a,K > the following inequalities are satisfied: 

\\u{t,-) - UooWs'^ ^ K — (^^^ Iho - ""oollsi- (57) 

Cl 

Proof. By Proposition 9.11 in [T21 Ch. 15], there is a unique solution u G C°°([0,oo) x S^). By the 
maximum principle, ||t((t, •)||5'i < HuqIIsi for t > 0. Define the function v = '^{u) of variables 
{t,x), where '^'{u) = k{u). The monotone function (p (of one variable) satisfies inequalities 

ci<ip' <C2. (58) 

The PDE ([56|)i reads as dtu = d'^^v, hence (in view of derivation dtv = ^'{u) dtu) is equivalent to 

dtv = k{u)d'^^v, f(0,x) = lp{uq{x)). 

Let us compare it with the linear PDE on S^, 

dtv = k{t,x) d'^^v, v{0,x) = ip{uo{x)), (59) 

where k{t,x) = k{u{t,x)) is given. Indeed, ci < k(t,x) < C2 for all t > and x £ S^. From the 
existence and uniqueness of a solution to (f59|) we conclude that v = v. Denote Uqo £ ^ the average 
of uq over 5"^, and set Voo = ^{uoo)- The function w = v{t,x) — Voo solves the linear PDE on 5^ 

dtw = k{t,x) d'^^w, w{0,x) = (f{uo{x)) — Voo, (60) 

By the theory of linear PDEs, ()60p possesses a fundamental solution G{t, x, y) which can be 
built by the classical parametrix method, and satisfies the inequalities < G < KG for some 
real K > 0. Here G is the fundamental solution of the heat equation dtu = a d'^^u for some 
constant a > 0, see survey [4]. Recall that G{t,x,y) = Ylj>o^~°''' Vj(^)0i(y)) where cpj denotes 
the eigenfunction (of operator —ad'^^ on 5^ with eigenvalue Xj = aj^) satisfying J^i (j)j{x) dx = 1. 
Indeed, {^{t, x) = J^i G{t, x, y) {ip{uo{y)) — Voo) dy. In particular, 

\\v{t,-) - < e~"*ll^^(0>-) - VooWs^- 

Thus, for all t > and x £ we have a priori estimate 

Mu{t,-)) -ip{uoo)\y < K e-'^^yiuo) - ip{uoo)\\s^- 
From this, using inequalities 

ci||M(t, •) - lioollsi < \\v{u{t, •)) - ip{Uoc)\\s^, e~°^\\ip{uo) - ip{Uoo)\\s^ < C2 e""*||iio - UocWs^, 
see §Si), we deduce ([57]). □ 
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Example 8. (i) Consider the heat equation over infinite (or semi-infinite) spatial interval M, 

dtu = d^^u + a{t,x), n(0, x) = no(x) G L^(R). (61) 
Here we assume that uq(x) is a bounded function. All the solutions of (|61|) . 

u{t,x)= / UQ{i)G{t,x,y)dy + / / a{y, s)G{t - s,x,y) dy ds, 



JR Jo JR 

include the heat kernel G{t,x,y) = (47rt)^"'^/^ e"^^"^-*^/^^*-* (the fundamental solution of homoge- 
neous (j6T|) for uq(x) = Sx{(,) - the Dirac delta function). For any function uq £ L2(M), a unique 
solution of the homogeneous heat equation, u{t,x) = J^uo{y)G{t,x,y) dy, converges uniformly to 
a linear function, as t — )• oo. Indeed, if uq is bounded then the linear function is constant. 

(ii) The rough Laplacian for tensors is defined by A = div V. For a unit circle S^, the eigenvalues 
of —A are the numbers \j = j"^ with eigenfunctions 4'j{x) = "^^6"'"'^) where j £ Z\ {0}. The heat 
equation on 5^, 

dtu = dl^u, u{0, •) = no G H^{S^), (62) 

has a unique solution for functions u e C([0, oo), H'^{S^)) fl C^((0,oo], L^(S'^)). By Sobolev em- 
bedding theorem, H^{S^) C G^{S^). The solution of ([62]) has the property that u{t,-) £ C°°(S'i) 
for all t > 0. Moreover, u{t, •) — )• no as t — >• oo, where uq = fgi uo{x) dx and 

\\u{t, •) — noil < e~*||no — no||. 

Consider the Jacobi theta function 

e(z, t) = = 1 + 2 y °° (e^^^)"' cos(27rnz), 

where = —1, z is a complex number and r is confined to the upper half-plane. Taking z = x gM. 
and r = Anti with real t > 0, we write 6{x,ATTti) = 1 + 2^^-^ e"^'^ " *cos(27rnx), which satisfies 
the heat equation (fU2|) i . Since lirn 9{x, Airt i) = J2nez ~ solution of can be specified 

by convolving the periodic boundary condition at t = with 9{x,4:TTti). 

(iii) Follow [3], denote U{t,x) = ^^^^i^^j^^it — 0) ^^"^ ^('") — T+u^- have a 3-parameter 
family u{t,x) = U{uj'^{t + /3i),a;(x + /32)) of exact solutions to the PDE ([56]) . dtu = dx{k{u)dxu)) 
on S^. Set = 1 and /3i = /32 = 0, hence u{t,x) = U{t,x) and no (2;) = ^cos^x+i ' -'^'^'^^^d' 
lim n(t, x) = lioo = for all x G S" . Since < n^ < 1, we conclude that ^ < k{u) < 1. Finally, we 
have ||n(t, •)||5i < e~* and UnoHgi = 1, that is consistent with ([ST]) . 

9.2 The generalized companion matrix 

Let Pn = — pik^^^ — ... — Pn-ik — Pn be a polynomial over R and A;i < ^2 < • • • < be the 
roots of Pn for n > 0. Hence, = (— l)*~^(Tj, where cij are elementary symmetric functions of the 
roots ki. The following generalized companion matrix (see [9j) plays a key role in this work: 



/ \ - \ 



I ••• 



- ^ 

n 

\(-l)"-lf<T„ ... -^a2 ai / 

Lemma 6 (see [9] ) . The matrix \63\l has the following properties: 

a) The characteristic polynomial of Bn,i is Pn- 

b) Vj = (1, 2kj, 3/c|, . . . ,n A;"~^) is the eigenvector of Bn,i for the eigenvalue kj. 

c) Bn^iV = VD, where V = { jA;*~"'^}i<jj<„ is the Vandermonde type matrix, and 

D = diag(A;i, . . . , kn) is a diagonal matrix. (If all ki 's are distinct, then V~^Bn^i V = D). 



(63) 
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Proposition 13 (see [9]). Let Ti(t,x) (i G N) be the power sums of smooth functions ki{t,x) (1 < 
i < n). Given m > 0, consider the infinite system of linear PDEs 

i 771 

^t'^i = TT i £ N. (64) 

2[i + m — 1) 



Then the n-truncated ^64\), i.e., Tn+i's are eliminated using suitable polynomials of ti, . . . ,Tn, is 



Remark that for m = 1 the system (j64p has diagonal form: dfTi = — ^ d^Ti, and for m = 2 the 
n-truncated system (fM|) reads as dtTi = — -Bn,i dxTi. 

Proposition 14. Let Ti{t,x) {i £ N) be as in Proposition \T3{. Given m> 0, consider the infinite 
system of linear PDEs 

dtTi = /"^ — - dl^n+rn^i, i e N. (65) 
2[i + m — 1) 

Then the n-truncated system l[U^) has the form 

Proof. By Proposition 1131 we have dxTn+i = Yl^=i^ij ^^Tj, where bij are the coefficients of the 
matrix Bn^m-i- Derivation of this, leads to d^xTn+i = Yl^=i ^ij '^"xx'^j + dxbij dxTj. □ 

Remark 5. For m = 1 the system (|65p has diagonal form: dtTi = ^d'^^Ti, and for m = 2 the 
n-truncated system ([65]) reads as 9(7" = Bn^i dxx^- 

For small values of m, m = 1, 2, the terms d'^^Tn+m are given by 

—y-r dlxTn+i = V" Gn-i+idlxTi + ai(T^, dxl^), (66) 

n + 1 ^-^1=1 I 

—^dlxTn+2 = {aian-i+1- (yn~-i+2)dlxTi + a2{'^ ,8x1^). (67) 

By Proposition I14t the last row of the matrix Bn^i or |(i?n,i)^ consists of the coefficients at d^^Ti's 
on the right-hand side of ([66|) . respectively, of (|67|) . and so on. 



Example 9. For fj = —2 5ji, (f65|) reduces to the (system of) heat equations 

dtTi = dxxTi, i G N, 

whose solution is known. Consider more complicated cases. 
1. For fj = —5j2, ([65]) reduces to the system 

dtn = dlxTi+i, ieN, (68) 

whose n-truncated version reads as: dtT^ = -Bn,i dxxT^ + 0'2{'^ , dx'r^)- For n = 2, we have two PDEs 

dtTi = ^ dlxT2, dtT2 = ^ dlxT3 = {t2 - Ti )<9^a;n + Tld^xT2, 

/ - \ 

and the matrix B21 = ( „ ^ I . If the roots ki 7^ k2, the eigenvectors of B2 1 are Vj = 

\ —2(T2 (Jl J 'J 

(1, 2 fcj), j = 1,2. For n = 3, (|68p reduces to the quasilinear system of three PDEs with the matrix 

i 

B^A =10 2 



>3,1 = I ^ ^ 3 

'2 



3(73 -|cJ2 CTl 
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whose eigenvalues are kj, and the eigenvectors are Vj 
2. For fj = —5j3, ([65]) reduces to the system 



3 i 

dtn = (i G N) ^ dt^ 

2(1 + 2) 



For n = 3, the system has the following matrix: 




3(T3 -^^2 O"! 

fo-lO-S i(o-3 - 0-1 0-2) - (T2) 







whose eigenvalues are |A;|. This series of examples can be continued as long as one desires. 
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